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Abstract 

The optimal patterns of synaptic conductances for spike generation in central neurons is a subject of con- 
siderable interest. Ideally, such conductance time courses should be extracted from membrane potential 
(V„,) activity, but this is difficult because the nonlinear contribution of conductances to the V„, renders 
their estimation from the membrane equation extremely sensitive. We outline here a solution to this prob- 
lem based on a discretization of the time axis. This procedure can extract the time course of excitatory 
and inhibitory conductances solely from the analysis of V,„ activity. We test this method by calculating 
spike-triggered averages of synaptic conductances using numerical simulations of the integrate-and-fire 
model subject to colored conductance noise. The procedure was also tested successfully in biological 
cortical neurons using conductance noise injected with dynamic-clamp. This method should allow the 
extraction of synaptic conductances from V^ recordings in vivo. 

Introduction 

Determining the optimal features of stimuli which aie needed to obtain a given response is of consider- 
able interest, for example in sensory physiology. Reverse-correlation is one of the most-used methods 
to obtain such estimates (Agiiera y Areas and Fairhall, 2003; Badel et al., 2006) and, in particular, the 
spike-triggered average (STA) is often used to determine optimal features linked to the genesis of action 
potentials (de Boer and Kuypers, 1968). The STA can be used to explore which feature of stimulus space 
the neuron is sensitive to, or to identify modes that contribute either to spiking or to the period of silence 
before the spike (Agiiera y Areas and Fairhall, 2003). Using intracellular recordings, it is straightforward 
to calculate the STA of the membrane potential (V,„), which yields the mean voltage trajectory preceding 
spikes. In contrast, it is much harder to determine the underlying synaptic conductance. Straightforward 
methods like recording at several different DC levels and estimating the total conductance from the ra- 
tio AI /AV fail, since the presence of a voltage threshold necessitates AV — >^ at the time of the spike, 
which, in turn, artificially suggests a divergence of the total conductance to infinity. Similarly, solving 
the membrane equation for excitatory and inhibitory conductances separately suffers from an additional 



complication: because the distance to threshold changes, the time courses of the average synaptic con- 
ductances depend on the injected current. 

Recent contributions (Badel et al., 2006; Paninski 2006a, 2006b) gave analytical expressions for the 
most likely voltage path, which in the low-noise limit approximates the STA, of the leaky integrate-and- 
fire (IF) neuron. In those cases, gaussian white noise current was considered as input. In Badel et al. 
(2006), a second state variable was added in order to obtain biophysically more realistic behavior. In 
Paninski (2006a, 2006b), in addition the exact voltage STA for the non-leaky IF neuron was computed, 
as well as the STA input current in discrete time. Here, a strong dependence of the STA shape on the 
time resolution dt was found without a stable limit as dt — > 0. It was argued heuristically that this 
behavior results from the fact that decreasing the time step corresponds to increasing the bandwidth of 
the input cun^ent, a point which was supported by numerical simulations (Paninski et al., 2004; Pillow 
and Simoncelli, 2003), in which a pre-filtering of the white noise input results in a stable limit STA. 

In this article, we focus on the problem of estimating the optimal conductance patterns required for 
spike initiation, based solely on the analysis of V„, activity. We consider neurons subject to conductance- 
based synaptic noise at both excitatory and inhibitory synapses. By discretizing the time axis, it is 
possible to obtain the probability distribution of conductance time courses that are compatible with the 
observed voltage STA. Due to the symmetry properties of the probability distribution, the STA time 
course of excitatory and inhibitory conductances can then be extracted by choosing the one with maxi- 
mum likelihood. We test this method in numerical simulations of the IF model, as well as in real cortical 
neurons using the dynamic-clamp technique, by comparing the estimated STA with the real STA deduced 
from the injected conductances. 

Material and Methods 

Models 

We considered neurons driven by synaptic noise described by two independent sources of colored con- 
ductance noise (point-conductance model (Destexhe et al., 2001)). The membrane equation of this sys- 
tem is given by: 
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Here, gi, ge{t) and gi{t) are the conductances of leak, excitatory and inhibitory currents, Vi, Vg, Vi are 
their respective reversal potentials, C is the capacitance and Idc ^ constant cuiTcnt. The subscript s in 
Eq. Q can take the values e, i, which in turn indicate the respective excitatory or inhibitory channel. We 
use gsQ and a, to indicate the mean and standard deviation (SD) of the conductance distributions, ^i (?) 
are gaussian white noise processes with zero mean and unit standard deviation. Throughout this article 
we use the correlation times = 2.728 ms and X; = 10.49 ms. 

This system was solved using numerical simulations of the leaky IF model, which was adjusted to 
match recordings of cortical neurons in slices (threshold —55 mV, refractory period 3 ms, reset —75 mV). 
Simulations were done using the NEURON simulation environment (Hines and Carnevale, 1997). To 
calculate STAs, approximately 1000 spikes occurring during spontaneous activity were used, each being 
preceded by a period of at least 100 ms of silence to avoid "contamination" of the V,„ STA by preceding 
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spikes. The same analysis protocols (see Results) were applied to the model and to experimental data. 

In order to address the influence of a dendritic filter on the reliability of our method, we used a two- 
compartment model based on that by Pinsky and Rinzel (1994). We removed all active channels and 
replaced them by an integrate-and-fire mechanism at the soma. The geometry (L = 3.18 /urn, diam = 
10 fjm) as well as the parameters not related to the spiking mechanism (gi = 10^^ S/cm^, Vl = —60 mV, 
capacitance c„, = 3 ^F/cm^, axial resistance = 5.87 x 10^ Hem) are identical for the two compart- 
ments. In addition, we chose a threshold for spiking V, = —55 mV at the soma. The parameters for 
leak conductance and capacitance needed for the estimation of the STAs of synaptic conductances from 
the V„,, g'^l' and C", were obtained by current pulse injection into the soma at the resting state. The 
superscript so indicates that these are effective values at the level of the soma. The values used were 
gf^ = 0.198 nS and = 5.86 pF. The parameters describing the distributions of synaptic conductances 
were chosen in a way such that the mean inhibitory conductance was four times that of excitation, and 
the latter was comparable to the leak conductance (geo = 0. 15 nS, gio = 0.6 nS). Standard deviations were 
assumed to be one third of the respective means (a^ = 0.05 nS, a, = 0.2 nS). 

In vitro experiments 

In vitro experiments were performed on 0.4 mm thick coronal or sagittal slices from the lateral portions of 
guinea-pig occipital cortex. Guinea-pigs, 4-12 weeks old (CPA, Olivet, France), were anesthetized with 
sodium pentobarbital (30 mg/kg). The slices were maintained in an interface style recording chamber at 
33-35 °C. Slices were prepared on a DSK microsUcer (Ted Pella Inc., Redding, CA) in a slice solution 
in which the NaCl was replaced with sucrose while maintaining an osmolarity of 307 mOsm. During 
recording, the slices were incubated in slice solution containing (in mM): NaCl, 124; KCl, 2.5; MgS04, 
1.2; NaHP04, 1.25; CaCb, 2; NaHCOs, 26; dextrose, 10, and aerated with 95% O2, 5% CO2 to a final pH 
of 7.4. Intracellular recordings following two hours of recovery were performed in deep layers (layer IV, 
V and VI) in electrophysiologically identified regular spiking and intrinsically bursting cells. Electrodes 
for intracellular recordings were made on a Sutter Instruments P-87 micropipette puller from medium- 
walled glass (WPI, IBFIOO) and beveled on a Sutter Instruments beveler (BV-IOM). Micropipettes were 
filled with 1.2 to 2 M potassium acetate and had resistances of 80-100 after beveling. 

The dynamic-clamp technique (Robinson et al., 1993; Sharp et al., 1993) was used to inject computer- 
generated conductances in real neurons. Dynamic-clamp experiments were run using the hybrid RT- 
NEURON environment (developed by G. Le Masson, Universite de Bordeaux), which is a modified ver- 
sion of NEURON (Hines and Camevale, 1997) running under the Windows 2000 operating system (Mi- 
crosoft Corp.). NEURON was augmented with the capacity of simulating neuronal models in real time, 
synchronized with the intracellular- recording. To achieve real-time simulations as well as data transfer 
to the PC for further analysis, we used a PCI DSP board (Innovative Integration, Simi Valley, USA) with 
4 analog/digital (inputs) and 4 digital/analog (outputs) 16 bits converters. The DSP boai^d constrains 
calculations of the models and data transfers to be made with a high priority level by the PC processor. 
The DSP board allows input (for instance the membrane potential of the real cell incorporated in the 
equations of the models) and output signals (the synaptic cuiTcnt to be injected into the cell) to be pro- 
cessed at regular- intervals (time resolution = 0.1 ms). A custom interface was used to connect the digital 
and analog inputs/outputs signals of the DSP board with the intracellular amplifier (Axoclamp 2B, Axon 
Instruments) and the data acquisition systems (PC -based acquisition software ELPHY, developed by G. 
Sadoc, CNRS Gif-sur-Yvette, ANVAR and Biologic). The dynamic-clamp protocol was used to insert 
the fluctuating conductances underlying synaptic noise in cortical neurons using the point-conductance 



model, similar to a previous study (Destexhe et al., 2001). According to Eq. ([T]) above, the injected cur- 
rent is determined from the fluctuating conductances and gi{t) as well as from the difference of the 
membrane voltage from the respective reversal potentials, loynciamp = —ge {V — K) — gi {V — Vi) ■ 

All research procedures concerning the experimental animals and their care adhered to the American 
Physiological Society's Guiding Principles in the Care and Use of Animals, to the European Council 
Directive 86/609/EEC and to European Treaties series no. 123, and was also approved by the local ethics 
committee "Ile-de-France Sud" (certificate no. 05-003). 



Results 

We first explain the method for extracting STAs from V,„ activity, then we present tests of this method 
using numerical simulations and intracellular recordings in dynamic-clamp. 



Method to extract conductance STA 

The procedure we follow here to estimate STA of conductances from V,„ activity is based on a discretiza- 
tion of the time axis. With this approach, a probability distribution can be constructed whose maximum 
gives the most likely conductance path compatible with the STA of the V,„. This maximum is determined 
by a system of linear equations which is solvable if the means and variances of conductances are known 
(for a method to estimate conductance mean and variance, see Rudolph et al., 2004). 

We start from the voltage STA, which is an average over an ensemble of event-triggered voltage 
traces. Its relation to the conductance STAs is determined by the ensemble average of Eqs. ([T]) and ([2]). 
In general, there is a strong correlation (or anti-correlation) between V{t) and gs{t) in time. However, 
it is safe to assume that there is no such correlation across the ensemble, since the noise processes ^s(f) 
corresponding to each realization are uncorrelated. Also, the ensemble average is commutative with the 
time derivative. Thus, we can rewrite Eqs. ([T]) and ^ to obtain 
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where = C/gi and (.),- denotes the ensemble average. In other words, the time evolution Eqs ([T]l and 
(|2]) also hold in terms of ensemble averages. In the following, we drop the bracket notation for legibility, 
but assume we are dealing with ensemble averaged quantities unless otherwise stated. 
We discretize Eq. ^ in time with a step-size At and solve for g^, 
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Since the series V'^ for the voltage STA is known, has become a function of ^g. In the same way, we 
solve Eq. (0]) for which have become gaussian distributed random numbers. 



There is a continuum of combinations {gg^^gf^' } that can advance the membrane potential from 
to V*^+^, each pair occumng with a probability 

T,- ( k+l k(- ' ^ 



where we have used Eq. Note that because of Eq. g'l and are not independent and is, 
thus, a unidimensional distribution only. Given initial conductances g^, we can now write down the 
probability p for certain series of conductances {gi}y=o,...,n to occur that reproduce a given voltage trace 

{V'}/=i,...,„+i: 

p = Y\pK (9) 

Due to the symmetry of the distribution p, the average paths of the conductances coincide with the most 
likely ones, so the cumbersome task of solving nested gaussian integrals can be circumvented. Instead, in 
order to determine the conductance series with extremal likelihood, we solve the n-dimensional system 
of linear equations 
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where X = Y!IZqX'^, for the vector {g^}. This is equivalent to solving = 0}k=i,...,n and involves 
the numerical inversion of an « x ?r-matrix. Since the system of equations is linear, if there is a solution 
for {^g}, plausibility arguments suggest that it is the most likely (rather than the least likely) excitatory 
conductance time course. The series {gf } is then obtained from Eq. 

Test of the accuracy of the method using numerical simulations 

To test this method, we first considered numerical simulations of the IF model in four different situations. 
We distinguished high-conductance states, where the total conductance is dominated by inhibition, from 
low-conductance states, where both synaptic conductances ai^e of compai^able magnitude. We also varied 
the standard deviations of the conductances such that for both high- and low-conductance states we 
have the cases a, > Og as well as > a,. The results are summarized in Fig. [H where the STA traces 
of excitatory and inhibitory conductances recorded from simulations are compared to the most likely 
(equivalent to the average) conductance traces obtained from solving Eq. (flOl ). In general, the plots 
demonstrate a very good agreement. 

To quantify our results, we investigated the effect of the statistics as well as of the broadness of the 
conductance distributions on the quality of the estimation. The latter is crucial, because the derivation 
of the most likely conductance time course allows for negative conductances, whereas in the simulations 
negative conductances lead to numerical instabilities, and conductances are bound to positive values. 
We thus expect an increasing error with increasing ratio SD/mean of the conductance distributions. We 
estimated the root-mean-square (RMS) of the difference between the recorded and the estimated con- 
ductance STAs. The results, summarized in Fig. |2l ai^e as expected. Increasing the number of spikes 
enhances the match between theory and simulation (Fig. |2]\ shows the RMS deviation for excitation, 

for inhibition) up to the point where the effect of negative conductances becomes dominant. In this 



example, where the ratio SD/mean was fixed at 0.1, the RMS deviation enters a plateau at about 7000 
spikes. The plateau values can also be recovered from the neighboring plots (i.e., the RMS deviations 
at SD/mean = 0.1 in Fig. |2t and D correspond to the plateau values in A and B). On the other hand, a 
broadening of the conductance distribution yields a higher deviation between simulation and estimation. 
However, at SD/mean = 0.5, the RMS deviation is still as low as ~ 2% of the mean conductance for 
excitation and ~ 4% for inhibition. 

To assess the effect of dendritic filtering on the reliability of the method, we used a two-compartment 
model based on that of Pinsky and Rinzel (1994), from which we removed all active channels and re- 
placed them by an integrate-and-fire mechanism at the soma. We repeatedly injected the same 100 s 
sample of fluctuating excitatory and inhibitory conductances in the dendritic compartment and performed 
two different recording protocols at the soma (Fig. [3K). First, we recorded in cuiTcnt-clamp in order to 
obtain the V,„ time course as well as the spike times. In this case, the leak conductance gf and the 
capacitance O''" were obtained from current pulse injection at rest. Second, we simulated an "ideal" 
voltage-clamp (no series resistance) at the soma using two different holding potentials (we chose the 
reversal potentials of excitation and inhibition, respectively). Then, from the currents /y^ and /v, , one can 
calculate the conductance time courses as 

where the superscript so indicates that these are the conductances seen at the soma (in the following 
refeiTcd to as somatic conductances). From these, we determined the parameters gj/^, g^Q, of and a™, 
the conductance means and standard deviations. In contrast to ge{t) and gi{t), the distributions of 
and gf'{t) are not Gaussian (not shown), and have lower means and variances. We compared the STA of 
the injected (dendritic) conductance, the STA obtained from the somatic V,„ using our method and the 
STA obtained using a somatic "ideal" voltage-clamp (see Fig.|3^-D), which demonstrated the following 
points: (1) as expected, due to dendritic attenuation, all somatic estimates were attenuated compared 
to the actual conductances injected in dendrites (compare light and dark gray curves, soma, with black 
curve, dendrite, in Fig. |3^-D); (2) the estimate obtained by applying the present method to the somatic 
Vm (dark gray curves in Fig.[3^-D) was very similar to that obtained using an "ideal" voltage-clamp at 
the soma (light gray curves). The difference close to the spike may be due to the non-Gaussian shape 
of the somatic conductance distributions, whose tails then become important; (3) despite attenuation, the 
qualitative shape of the conductance STA was preserved. We conclude that the STA estimate from V„t 
activity captures rather well the conductances as seen by the spiking mechanism. 

Test of the method in real neurons 

We also tested the method on voltage STAs obtained from dynamic-clamp recordings of guinea-pig 
cortical neurons in slices. In real neurons, a problem is the strong influence of spike-related voltage- 
dependent (presumably sodium) conductances on the voltage time course. Since we maximize the global 
probability of ge{t) and gi{t), the voltage in the vicinity of the spike has an influence on the estimated 
conductances at all times. As a consequence, without removing the effect of sodium, the estimation 
fails (see Fig. |4l). Fortunately, it is rather simple to con^ect for this effect by excluding the last 1-2 ms 
before the spike from the analysis. The corrected comparison between the recorded and the estimated 
conductance traces is shown in Fig. [5] 

Finally, to check for the applicability of this method to in vivo recordings, we assessed the sensitivity 



of the estimates with respect to the different parameters by varying the values describing passive prop- 
erties and synaptic activity. We assume that the total conductance can be constrained by input resistance 
measurements, and that time constants of the synaptic currents can be estimated by power spectral analy- 
ses (Destexhe and Rudolph, 2004). This leaves gi, C, geo, Oe and Oi as the main parameters. The impact 
of these parameters on STA conductance estimates is shown in Fig. [6l Varying these parameters within 
±50% of their nominal value led to various degrees of error in the STA estimates. The dominant effect 
of a variation in the mean conductances is a shift in the estimated STAs, whereas a variation in the SDs 
changes the curvature just before the spike. 

To address this point further, we fitted the estimated conductance STAs with an exponential function: 

Ut)=G,{l+K,e'-^), (12) 

where s again takes the values e, i for excitation and inhibition, respectively, fo is chosen to be the time at 
which the analysis stops. Fig.|7]gives an overview of the dependence of the fitting parameters Ge, Gj, T^. 
and Ti on the relative change of gi, geo, Oe, and C. For example, a variation of geo has a strong effect 
on Ge and G,-, but affects to a lesser extent Tg and 7]-, while the opposite was seen when varying and 

Discussion 

Understanding the transfer function of a neuron from synaptic input to spike output would ideally require 
the simultaneous monitoring of both the synaptic conductances and the cell's firing. CuiTcnt methods for 
extracting synaptic conductances rely on intracellular- recordings performed at different holding poten- 
tials (in voltage-clamp) or different current levels (in current-clamp; e.g. Borg-Graham et al. 1998) and, 
as a consequence, they do not allow the establishment of a direct correspondence between synaptic con- 
ductances and spikes. Although these methods have been very useful, for example in establishing the 
synaptic structure of sensory receptive fields in a variety of systems (Monier et al., 2003; Wehr^ and 
Zador, 2003; Wilent and Contreras, 2005), they do not distinguish between trials that effectively produce 
spikes at a given latency and those that do not. 

Here, we have presented a method to extract the average excitatory and inhibitory conductance pat- 
terns directly related to spike initiation. As illustrated in Fig. [H this method can extract spike-related 
conductances based solely on the knowledge of V,„ activity. First, the STA of the V,„ is computed from 
the intracellular recordings. Next, by discretizing the time axis, one estimates the "most likely" con- 
ductance time courses that are compatible with the observed STA of V„,. Due to the symmetry of their 
distribution, the average conductance time courses coincide with the most likely ones, so integration 
over the entire stimulus space (whose dimension depends on the STA interval as well as on the temporal 
resolution) can be replaced by a differentiation and subsequent solution of a system of linear equations. 
Solving this system gives an estimate of the average conductance time courses. We demonstrated that 
this estimation gives reasonably accurate estimates for the leaky IF model, as well as in real neurons 
under dynamic-clamp. 

Like any other method, this method suffers from several sources of error. Errors can result from non- 
linearities in the I-V curve of the neuron, e.g. those due to voltage-dependent conductances. In agreement 
with this, we have shown that the subthreshold activation of spike-generating currents close to thr^eshold 
can lead to severe misestimations of the conductances (Fig. |4ll. This problem can be circumvented by 
excluding a short period (1-2 ms) preceding the spike. To avoid contamination by voltage-dependent 
currents, this method should be complemented by a check for I-V curve linearity in the range of V„j con- 



sidered. Note that a linear I-V curve does not guarantee the absence of voltage-dependent conductances. 
For example, if the mean interspike interval of the cell becomes too short, spike-related potassium cur- 
rents might be present during a substantial fraction of the STA interval and could affect the estimation. 
This might diminish the applicability of the method to neurons spiking at high frequency, in particular 
to fast-spiking intemeurons. Also, strong subthreshold dendritic conductances that are very remote from 
the soma could influence the STA estimate without being visible in the I-V curve. On the other hand, in 
cases where it is possible to parameterize these nonlinearities, they can be included in Eq. |5] It should 
thus be possible to extend the method in order to apply it to more complex models, for example the ex- 
ponential integrate-and-fire model (Fourcaud-Trocme et al., 2003). Another possible extension would be 
to include voltage-dependent terms such as N-methyl-D-aspartate (NMDA) receptor-mediated synaptic 
currents, although such currents probably have a limited contribution at the range of considered here 
(below -50 mV). 

Another source of error may arise from "negative conductances". The present model of synaptic 
noise assumes that conductances are Gaussian-distributed, but if the standard deviation becomes com- 
parable to the mean value of the conductances, the Gaussian distribution will include negative conduc- 
tances, which are unrealistic. This is an important limitation of representing synaptic conductances by 
Gaussian-distributed noise ("diffusion approximation"). However, this type of approximation seems to 
apply well to cortical neurons in vivo, which receive a large number of inputs (Destexhe et al., 2001). 
In vivo measurements so far indicate that the standard deviation is much smaller than the mean for both 
excitatory and inhibitory conductances (Haider et al., 2006; Rudolph et al., 2005), which also indicates 
that the diffusion approximation is valid in this case. Such a check for consistency is a prerequisite for 
applying the present method. 

Previous work related to the question of spike-tiiggered stimuli was mainly focused on white noise 
current inputs, and showed that no stable finite input average exists in the limit dt —^0 (Paninski 2006a). 
Other contributions shed light on the question of membrane voltage STAs, for the leaky IF neuron as well 
as for biophysically more plausible models. However, so far no procedure was proposed to solve this 
problem of reverse correlation for conductance noise inputs. The method we propose here attempts to fill 
this gap, and directly provides a procedure that can be applied to real neurons. To this end, the present 
method must be complemented by measurements of the mean and standai^d deviation of excitatory and 
inhibitory conductances. Such measurements can be obtained either by voltage-clamp (Haider et al., 
2006), or by current-clamp, as recently proposed (Rudolph et al, 2004, 2005). Combining the latter 
method with the present method, it should now be possible to directly extract conductance patterns from 
Vm recordings in vivo, and thus obtain estimates of the conductance variations related to spikes during 
natural network states. 
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Figure 1 : Test of the STA analysis metliod using an IF neuron model subject to colored conductance noise. A. 
Scheme of the procedure used. An IF model with synaptic noise was simulated numerically (bottom) and the 
procedure to estimate STA was applied to the V,„ activity (top). The estimated conductance STAs from V„, were 
then compared to the actual conductance STA in this model. Bottom panels; STA analysis for different conditions, 
low-conductance states (B,C), high-conductance states (D,E), with fluctuations dominated by inhibition (B,D) 
or by excitation (C,E). For each panel, the upper graph shows the voltage STA, the middle graph the STA of 
excitatory conductance, and the lower graph the STA of inhibitory conductance. Solid lines (grey) show the average 
conductance recorded from the simulation, while the dashed line (black) represents the conductance estimated from 
the V„,. Parameters in B: geQ=6 nS, gio=6 nS, Of,=0.5 nS, 0,=1.5 nS; C: geo=6 nS, gio=6 nS, 0£.= 1.5 nS, a,=0.5 nS; 
D: geo=20 nS, §,o=60 nS, 0^=4 nS, o,=12 nS; E: geo=20 nS, gm=60 nS, 0^=6 nS, o,=3 nS 
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Figure 2: The root-mean-square (RMS) of the deviation of the estimated from the recorded STAs. A: RMS 
deviation as a function of the number of spikes for the STA of excitatory conductance, where the SD of the 
conductance distribution was 10% of its mean. The RMS deviation first decreases with the number of spikes, but 
saturates at ~ 7000 spikes. This is due to the effect of negative conductances, which are excluded in the simulation 
(cf. C). B: Same as A for inhibition. C: RMS deviation for excitation as a function of the ratio SD/mean of the 
conductance distribution. The higher the probability of negative conductances, the higher the discrepancy between 
theory and simulation. However, at SD/mean = 0.5, the mean deviation is as low as ^ 2% of the mean conductance 
for excitation and 4% for inhibition. D: Same as C for inhibition. 
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Figure 3: Test of the method using dendritic conductances. A. Simulation scheme: A 100 s sample of excitatory 
and inhibitory (frozen) conductance noise was injected into the dendrite of a 2-compartment model (1). Then, two 
different recording protocols were performed at the soma. First, the Vm time course was recorded in current-clamp 
(2), second, the currents corresponding to two different holding potentials were recorded in voltage-clamp (3). 
From the latter, the excitatory and inhibitory conductance time courses were extracted using Eq.[TT] B. STA of total 
conductance inserted at the dendrite (black), compared with the estimate obtained in voltage-clamp (light gray) 
and with that obtained from somatic V,„ activity using the method (dark gray). Due to dendritic attenuation, the 
total conductance values measured are lower than the inserted ones, but the variations of conductances preceding 
the spike are conserved. C. Same as B, for excitatory conductance. D. Same as B, for inhibitory conductance. 
Parameters: 0.15 nS, gm = 0.6 nS, a, = 0.05 nS, a,- = 0.2 nS, =0.113 nS, gj^^ = 0.45 nS, of = 0.034 
nS, = 0.12 nS, where the superscript so denotes quantities as seen at the soma. 
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Figure 4: The effect of the presence of additional voltage-dependent conductances on the estimation of the 
synaptic conductances. Gray, solid lines indicate recorded conductances; black, dotted lines indicate estimated 
conductances. In this case, the estimation fails. The sharp rise of the voltage in the last ms before the spike 
requires very fast changes in the synaptic conductances, which introduces a considerable error in the analysis. 
Parameters used: geQ=32 nS, gio=96 nS, 0£.=8 nS, 0,=24 nS. 
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Figure 5: Test of the method in real neurons using dynamic-clamp in guinea-pig visual cortical slices. A. Scheme 
of the procedure. Computer-generated synaptic noise was injected in the recorded neuron under dynamic-clamp 
(bottom). The V,„ activity obtained (top) was then used to extract the STA of conductances, which was compared 
to the STA directly obtained from the injected conductances. B. Results of this analysis in a representative neuron. 
Black lines show the estimated STA of conductances from V,„ activity, grey lines show the STA of conductances 
that were actually injected into the neuron. The analysis was made by excluding the data from the 1.2 ms before 
the spike to avoid contamination by voltage-dependent conductances. Parameters for conductance noise were as 
in Fig. H 
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Figure 6: Deviation in the estimated conductance STAs in real neurons using dynamic-clamp due to variations 
in the parameters. The black lines represent the conductance STA estimates using the correct parameters, the 
gray areas are bound by the estimates that result from variation of a single parameter (indicated on the right) by 
± 50 %. Light gray areas represent inhibition, dark gray areas represent excitation. The total conductance (leak 
plus synaptic conductances) was assumed to be fixed. A variation in the mean values of the conductances evokes 
mostly a shift in the estimate, while a variation in the SDs influences the curvature just before the spike. 
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Figure 7: Detailed evaluation of the sensitivity to parameters. The conductance STAs were fitted with an expo- 
nential function fs{t) = Gs{l +Ksexp{{t — t(j)/Ts), s — e,i. to is chosen to be the time at which the analysis stops. 
Each plot shows the estimated value of Ge, Gi, Tg or Ti from this experiment, each curve represents the variation 
of a single parameter (see legend). 
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Figure 8: Scheme of the method to extract spike-triggered average conductances from membrane potential 
activity. Starting from an intracellular recording (top), the spike-triggered average (STA) membrane potential (V,„) 
is computed (leftmost panel). From the STA of the V,„, by discretizing the time axis, it is possible to estimate the 
STA of conductances (bottom) by maximizing a probability distribution (see text). This step requires knowledge 
of the values of the average conductances and their standard deviations (geo, gio, Oe, O,, respectively), which must 
be extracted independently (rightmost panel). 



